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Abstract 

Specific heat C(T) measurements were made on single crystals of the superconducting filled 
skutterudite series Pr(Osi_a;Ru x )4Sbi2 down to 0.6 K. Crystalline electric field fits in the normal 
state produced parameters which were in agreement with previous measurements. Bulk super- 
conductivity was observed for all values of the Ru concentration x with transition temperatures 
consistent with previous experiments, confirming a minimum in T c at x = 0.6. The C(T) data 
below T c appear to be more consistent with power law behavior for x = (PrOs4Sbi2), and with 
exponential behavior for 0.05 < x < 0.2. An enhanced electronic specific heat coefficient 7 was 
observed for x < 0.4, further supporting x ~ 0.6 as a critical concentration where the physical 
properties abruptly change. Significant enhancement of AC/T C above the weak coupling value was 
only observed for x = and x = 0.05. 

PACS numbers: 65.40.Ba, 71.27.+a, 74.25.Bt, 74.62.Dh 
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I. INTRODUCTION 



The filled skutterudite compound PrOs4Sbi2 has proven to be an intriguing and unusual 
material since its discovery as the first Pr-based heavy fermion superconductor a few years 
ago.i*2 Both the heavy fermion (m* ~ 50 m e ) and the superconducting (T c = 1.85 K) states 
display very unusual properties. The ground state of the Pr 3+ ion in PrOs4Sbi2 that arises 
from the splitting of the Pr 3+ J = 4 multiplet in a crystalline electric field (CEF) is nonmag- 
netic, and is either a Ti singlet or a T3 doublet. The Ts triplet first excited state is ~ 10 K 
above the ground state, with the other excited states following at ~ 100 K (r 4 triplet) and 
~ 300 K (r 3 or r 1; respectively). In the superconducting state, PrOs 4 Sbi 2 exhibits multiple 
transitions in specific heat^* 3 - and magnetic penetration depth,— and may also contain multi- 
ple superconducting phases.^ The nature of the superconducting energy gap is also not clear: 
muon spin rotation^ (/zSR) and Sb-nuclear quadrupole resonance^ (Sb-NQR) measurements 
indicate isotropic and strong-coupling superconductivity, tunneling spectroscopy measure- 
ments support a nearly fully gapped but unconventional superconducting order parameter,— 
and data from thermal conductivity in a magnetic field^ and magnetic penetration depth 9 
are consistent with point nodes in the energy gap. Additional /iSR measurements reveal 
possible time-reversal symmetry breaking in the superconducting state, further suggesting 
that the superconducting state does not have s-wave symmetry.— 

PrRu 4 Sbi2 is a much simpler compound than PrOs 4 Sbi 2 . It is also superconducting 
(T c = 1.1 K), but displays more conventional properties. ^ From Sb-NQR measurements, 
the superconductivity appears to be weak-coupling with an isotropic energy gap. 12 Mag- 
netic penetration depth measurements yield moderate coupling and a fully gapped order 
parameter.— In addition, PrRu 4 Sbi2 is not a heavy fermion compound; it has an electronic 
specific heat coefficient 7 ~ 10 times smaller than that of PrOs 4 Sbi2- Features in the physi- 
cal properties of PrRu 4 Sbi 2 could be described by a CEF model with a T 1 ground state and 
a T 4 first excited state separated by ~ 70 K . 11 ' 14 

The Pr(Osi_ 3; Ru a: ) 4 Sbi2 series of compounds was previously studied through measure- 
ments of magnetic susceptibility x(T) and electrical resistivity p(T).— Superconductivity 
was found to persist for all values of the Ru concentration x, with a minimum in the T c — x 
phase diagram at x = 0.6 where T c = 0.75 K. This minimum may arise from a competition 
between the heavy fermion superconductivity of PrOs 4 Sb 12 (x = 0) and the BCS supercon- 
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ductivity of PrRii4Sbi2 (x = 1). Based on theoretical models, it has recently been suggested 
that there may be a mixed-parity superconducting state near this minimum in T c }- CEF 
effects were also observed for all values of x in the normal state of x(T) and p(T), with 
the splitting between the ground state and the first excited state increasing monotonically 
with x between x = and x = 1. For x(T), fits with a T 3 ground state were consistent 
with the data for all values of x, while fits with a I\ ground state were only satisfactory 
near the extremal values of the Pr(Osi_ x Ru x )4Sbi2 series. The p(T) data were also fit with 
CEF equations, and although the fits were insensitive to the degeneracy of the ground state, 
they were still able to provide level splittings consistent with those derived from the x{T) 
data. In the present study, the specific heat C(T) of P^Osi^Ru^Sb^ was measured down 
to ~ 0.6 K, to further investigate the normal and superconducting state properties of this 
extraordinary system. 

II. EXPERIMENTAL DETAILS 

The single crystal specimens of Pr(Osi_ x Ru x ) 4 Sbi 2 investigated in this work were identi- 
cal to those previously studied.— Specific heat C was measured as a function of temperature 
between 0.6 and 50 K in a 3 He semiadiabatic calorimeter by using a standard heat pulse 
technique. The samples were attached to a sapphire platform with Apiezon N grease. The 
data presented in this work were taken from experiments on collections of single crystals with 
total masses between 11 and 114 mg. X-ray measurements show no signs of multiple phases 
in the doped materials; however, there was some sample dependence of the superconducting 
transition in electrical resistivity, especially on the doped materials.— Thus, the C(T) data, 
especially at the superconducting transitions, are expected to be slightly broadened by the 
sample dependence of the crystals. 

III. RESULTS AND DISCUSSION 

Displayed in Fig. [Hare specific heat divided by temperature C/T vs T data for various Ru 
concentrations x of the Pr(Osi_ x Ru x ) 4 Sbi 2 series for temperatures between 0.6 and 20 K. 
The maximum of the Schottky anomaly, associated with the CEF splitting of the Pr 3+ energy 
levels, noticeably decreases in magnitude and shifts to higher temperatures with increasing 
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x. All of the C{T) data for Pr(0si_ x Ru x )4Sbi2 were fitted between their superconducting 
transitions T c and 10 K by an equation including electronic, lattice, and Schottky terms: 



(where 0£> is the Debye temperature), and C^iT) is the Schottky specific heat anomaly for 
a two level system arising from the energy difference between the CEF ground state and the 
first excited state, scaled by a factor r. The results of these fits are listed in Table These 
fits find values of 0£> for the end member compounds comparable to other single crystal 
results of 165 K for PrOs 4 Sbi 2 17 and 232 K for PrRu 4 Sbi 2 .^ 

The Schottky specific heat anomaly Csch(^) for a two- level system is given by 



where 5 is the energy difference in units of K between the two levels, and go and g\ are the 
degeneracies of the ground state and excited state, respectively— In zero magnetic field, this 
equation is independent of whether or not the local symmmetry of the Pr 3+ ions is cubic 
or tetrahedral. It was found during the fitting procedure that in order for the fits to be 
accurate, one or more of the terms in Eq. ^ had to be scaled. Fits were done both with 
r modifying the Schottky term, and with r multiplying the entire equation. The former 
case would be interpreted as some internal broadening of the energy levels or as an overall 
transfer of entropy to the itinerant electrons due to hybridization, while the latter case would 
imply impurity phases (most likely free Sb) causing an overall overestimate of the sample 
mass. While both possibilities produced good qualitative fits, the values for 7 resulting from 
assuming an overall scaling were extremely large and not physically reasonable. Therefore, 
all the fits presented here were exactly as shown in Eq. Q with r only modifying the Schottky 
anomaly term. 

The normal state fits were only performed up to 10 K so that the Ciattice ~ (3T 3 approx- 
imation would more likely be accurate; however, the lattice terms are clearly the smallest 
in this temperature range compared to the other terms, and are thus difficult to accurately 
fit. This appears to especially be true for x = 0.2 and x = 0.4, where Qd is suppressed 
compared to the end member compounds, and which may be due to the disorder inherent 
in the substituted compounds. Because of the uncertainty in the accuracy of the fit values 



C/T = 1 + f3T 2 + rC Sc ^/T. 



(1) 



Here 7 is the electronic specific heat coefficient, (3 cx D 3 is the lattice specific heat coefficient 
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for 0£>, the error in the other parameters was estimated by varying Qd by ±10 and refitting 
the data. The largest error in the normal state arises in 7, while the errors in 5 and r are 
much smaller. These errors are represented in Table and in the figures where appropriate. 

The question of whether or not the ground state in PrOs4Sb 12 is a T 3 doublet or a Ti 
singlet has been contentious since the heavy fermion superconductivity of PrOs4Sbi2 was 
discovered. In our original reports of heavy fermion superconductivity in PrC^Sbi^r^ fits to 
the magnetic susceptibility yielded two possible Pr 3+ crystalline electric field splittings, both 
with a T 5 first excited state and either a Ti singlet ground state or a T3 nonmagnetic doublet 
ground state with a quadrupole moment. At the present time, it appears that the overall 
data are better explained by a Ti singlet ground state . 19 i 20 i 21 Nevertheless, it was felt that fits 
to the C(T) data for both possibilities should be made. As can be seen from Table and Fig. 
|2[b), when fitting the normal state data up to 10 K, both possible ground state fits result 
in reasonable values that agree well with those previously published for Pr(Osi_ x Rua;)4Sbi2 
based on x(T) an d p(T) measurements.— However, for x < 0.4, the scaling factor r is much 
closer to 1 for a T3 ground state compared to Y\. In fact, the best T3 ground state fit for 
x = 0, pure PrOs4Sbi2, results in almost no scaling whatsoever. (The present fits were 
performed on data from a measurement of single crystals, while some previously published 
fits were performed on data from a measurement of a pressed pellet . 1 ^ 22 ) These results of 
r-p 3 = 1.01 and rp 1 = 0.56 for x = are very consistent with fits performed on data from 
measurements of single crystals by Vollmer et al, which resulted in the values rr 3 ~ 0.99 
and rpj ~ 0.5.^ As x increases, the scaling factor decreases, indicating a suppression of 
the Schottky anomaly. At face value, the fact that r-p 3 is always closer to 1 than r-p 1 for 
x < 0.4 could be considered as support for a T 3 ground state. It has been suggested that 
the suppression of the Schottky anomaly, especially for a Ti ground state fit for PrOs4Sbi2, 
could result from an energy dispersion due to Pr-Pr interactions^ or hybridization between 
the Pr f-electrons and ligand states.— As these arguments can apply equally well to either 
ground state, these normal state results for x < 0.4 appear unable to discern between Ti 
and T3 ground states. 

For x > 0.4, the scaling factor r increases rapidly. PrRu4Sb 12 exhibits a situation com- 
plementary to that of PrOs4Sbi2: while rp 3 is still 1.8 — 1.9 times larger than rp i; it is r-p t 
that is closer to 1. The accuracy of these results could be affected by the temperature limits 
of the fits; for these large splittings, the maximum of the Schottky anomaly is well above 10 
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K. However, the calculated values of S for the different ground states agree very well overall 
with the values previously measured,— supporting the monotonic increase of the splitting 
between the ground state and the first excited state throughout the Pr(Osi_a;Ru x )4Sbi2 
series (Fig. I2fb)). 

In order to get a more accurate determination of T c and the specific heat jump at T c , 
AC/T C , the CEF and lattice fit results were subtracted from the data at low temperatures, 
leaving only the electronic specific heat. The subtractions were carried out using both the 
T 3 and the I\ ground state CEF fits, including those due to the variation of 6^ in order 
to further estimate the error of the parameters. All the subtractions resulted in exactly 
the same T c values, and very similar (within experimental error) AC/T C values, which are 
both listed in Table ITT1 These values of T c are plotted as a function of x in Fig. |2fa), along 
with previously measured values; the error bars for the T c points represent the width of the 
transitions in the respective measurements. All the data agree very well, and the minimum 
at x = 0.6 is also reproduced in the current data. With the exception of pure PrOs4Sbi2 
(i.e., for x > 0), the shapes of the measured C(T) curves below T c were nearly the same for 
either Ti or r 3 -based subtraction. However, for PrOs4Sbi2, a significant difference can be 
seen in the slope of the data below T c for the two different subtractions. The C(T) data 
after subtraction of Ci at (T) and Csch(^) for a Ti or T 3 ground state, C e \(T), are shown in 
Fig. 01 All the concentrations are shown for the Ti ground state subtraction, and the T 3 
ground state subtraction is also shown for x = 0. 

The data below T c , after the lattice and CEF terms were subtracted, were fit to both 
power-law and exponential functions, for energy gaps with and without nodes, respectively. 
These functions are typically considered to be valid only at very low temperatures. However, 
there are several examples of heavy fermion superconductors which appear to display power- 
law behavior up to near T c (e.g. , 25 i 26 ). The current experiment had a low-temperature limit of 
0.6 K, which is effectively a base temperature due to the large nuclear Schottky contribution 
at lower temperatures . 1 ^ 2 ? The high temperature limit of the fit was chosen to be |T C in 
light of the above referenced examples, and also to avoid possible spurious effects due to the 
width of the superconducting transitions. Because of these constraints, only the samples 
with < x < 0.2 could be fitted below T c , as there were not enough data points at the other 
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concentrations to give a reliable fit. The power- law fit was of the form 



C/T = 7& p + AT n , 



(3) 



suggested for energy gaps containing nodes, 27 and the exponential fit was of the BCS form^ 



where 7 S represents an electronic specific heat coefficient in the superconducting state, A 
and B are fitting constants, and A e is proportional to the energy gap in the BCS theory of 
superconductivity. The results of the application of these fits to the x = and x = 0.05 data, 
for both ground state subtractions, are shown in Fig. where the solid lines correspond 
to the power law fits and the dashed lines represent the exponential fits. For x — 0, the 
power-law fits extrapolate from the highest fit temperature of ~ 1.2 K all the way up to the 
transition, while the exponential fits deviate from the data right above 1.2 K. For x = 0.05 
and higher, the converse is true, as the exponential fits extrapolate to the higher temperature 
data (above ~ 1.1 K for x = 0.05 and x = 0.1, and above ~ 1.0 K for x = 0.2) much 
more accurately than the power law fits. While these extrapolations cannot be taken by 
themselves as proof of the superiority of one fit over the other, they are certainly suggestive 
and intriguing. The broadened superconducting transitions in the Ru-doped samples in 
particular may make the exponential fit appear more appropriate. 

The constants from the fits below T c are also listed in Table |H] The listed errors are due 
to the variation of Qd in the normal state fits. The error in T c was taken to be the width of 
the superconducting transition and is presented graphically in Fig. |2^; the error in AC/T C 
from the variation of was negligibly small compared to the error inherent in the equal 
area construction for determining AC/T C , which is represented graphically in Fig. 03 Since 
these fits are phenomenological, the absolute values of the resulting parameters should not 
necessarily be trusted. However, comparing the fits and the samples to one another can 
prove instructive. For pure PrOs4Sbi2, the ground state subtraction results in fits that 
are much different from the Ti ground state subtraction results. The other concentrations 
have similar fit parameters for the two different ground states. In addition, the values for 
7f, n, and A e , and the errors associated with them, are much larger for the r% x = fits 
than for the other T 3 fits, falling well outside the spread of the other three data points. If 
the fits are indeed accurate, then the value of n for the Ti x = data is comparable with the 




(4) 
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n = 2 expected in C/T for a heavy fermion compound with point nodes in the energy gap.— 
The values of A e for the T\ data of A e /T c « 2.2 — 2.6 are moderately enhanced compared 
to the weak-coupling value of A e /T c = 1.76 (for A e in units of K). The parameter 7 S can be 
interpreted as a portion of the sample that is normal or gapless; however, the values of 7 S 
could also simply be artifacts of the fit, especially in the case of the power law fit for x = 0.1. 
The large discrepancy between the fits below T c for PrOs4Sbi2, even if the absolute values 
are not entirely accurate, do lend phenomenological support for T\ being the ground state 
in this compound. It is interesting that the analysis of C(T) in the superconducting state 
for samples containing a small amount of Ru shows only small sensitivity to the choice of 
Ti or T 3 for the Pr 3+ ground state. 

Fig. shows a comparison between the electronic specific heat coefficient 7 calculated 
from the fits of the normal state C/T data and a value estimated from AC/T C . In the BCS 
theory of superconductivity, AC /^T c = 1.43. This is not likely to be true in the case of the 
unconventional superconductivity in PrOs4Sb 12 . However, it is expected that there will be 
some proportional relation between AC/T C and 7 (e.g., AC/T C oc 7), and so it can still be 
instructive to view this graphically. It can be seen in Table HTl and Fig. that 7(2;) derived 
from the normal state fits exhibits a peak at x — 0.05, decreases with x to a minimum at 
x = 0.6, and then slowly rises with x to x — 1. In contrast, AC/T C starts out extremely 
large for x = and decreases very quickly, again to a minimum value at x = 0.6. The 
"hump" in the data at x = 0.4 is due to the extremely broad superconducting transition at 
this concentration. 

The discrepancy between the values of 7 determined from the normal state data and 
calculated from the superconducting transition is an interesting one. The minimum in 7 at 
x = 0.6 for the measured 7, along with the minimum in T c , strongly suggests that something 
unusual is happening with the physical properties at this concentration. In contrast, the 
enhanced AC/T C for x = and x = 0.05 imply that strong-coupling superconductivity is 
only present for these two concentrations. Taken at face value, this could mean that for 
0.1 < x < 0.4, the heavy electrons are not participating in the superconductivity, which 
would imply that the superconductivity is nearly conventional for x > 0.1 (and possibly 
also for x = 0.05 if the exponential behavior below T c is appropriate). A more conventional 
explanation could be that the Ru substitution broadens the superconducting transitions 
in C(T) enough to result in an underestimate for AC/T C . This would only appear to be 
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the case for 0.05 < x < 0.4, however, as the superconducting transitions for x = 0.6 and 
x = 0.85 are as sharp or sharper than for x — 1. It would be unusual for disorder to play a 
strong role only for low Ru concentrations. 

In Figs. El and 0] it can be seen that the sharp double superconducting transition present 
in pure PrOs 4 Sbi 2 is not obviously apparent in the Ru-doped samples, although this may 
be obscured by the width of the transitions. Chia et al. have recently suggested that this 
double transition may arise from two superconducting phases which react differently to Ru 
substitution; i.e., one is unconventional and is quickly destroyed by impurities, while the 
other is more conventional and persists throughout the entire series.— It will be interesting 
to see whether or not this conjecture can tie together the seemingly contradictory results 
obtained through experiments over the past few years. Indeed, in studies on the specific 
heat of La-doped PrOs4Sbi2, preliminary results suggest that only the lower transition of 
PrOs 4 Sbi2 is suppressed for small La concentrations.— This could also be the case for the 
present Ru substitution studies. On the other hand, measurements of magnetic penetration 
depth suggest two-band superconductivity in PrOs4Sbi2, where the two bands are coupled 
by Josephson pair tunnelling, which requires the two order parameters to have the same 
symmetry- The difference in the behavior below T c between PrOs4Sb 12 and the Ru-doped 
samples could perhaps be accounted for by the suppression of one of the superconducting 
bands. Measurements on samples with small Ru concentrations are underway in order to 
more closely track the evolution of the superconducting properties in this remarkable series 
of compounds. 



IV. SUMMARY 



The specific heat of single crystal samples of Pr(Osi_ x Ru x )4Sbi2 was measured down to 
0.6 K. Fits to the normal state data resulted in CEF parameters that monotonically increase 
with x, in agreement with previously reported data. The superconducting transitions in 
C(T) were also consistent with previously reported x(T) and p(T) data, confirming the 
minimum in T c at x = 0.6. Below T c , a power law fit possibly corresponding to point nodes in 
the energy gap could better describe the data for PrOs4Sbi2, but an exponential fit associated 
with an isotropic energy gap was more appropriate for the data with 0.05 < x < 0.2. The 
electronic specific heat 7 was inferred from both normal state data and the ratio AC/T C . The 
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normal state fits revealed an enhanced 7 for x < 0.4, reaching a minimum value at x = 0.6, 
the same concentration as the minimum in T c . However, AC/T C was only significantly 
enhanced for x = and x = 0.05. 
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TABLE I: Physical properties of samples of P^Osi-^Ru^Sb^, determined from normal state 
specific heat data. The parameter 5 is the splitting between the ground state and the first excited 
state in the Schottky anomaly, r is the scaling factor for the Schottky anomaly, 7 is the estimated 
electronic specific heat coefficient, and Qd is the estimated Debye temperature. The errors in the 
parameters were determined by allowing 6d to vary by ±10 K within the fits (see text for details). 
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FIG. 1: Specific heat divided by temperature C/T below 20 K for single crystal samples of 
Pr(Osi_ z Ru :c )4Sbi2. 
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TABLE II: Physical properties of samples of Pr(Osi_ I Ru a; )4Sbi2, determined from superconducting 
state specific heat data. The parameter T c is the superconducting transition temperature, AC is 
the jump in C(T) at T c , 7 S is the electronic specific heat coefficient in the superconducting state, 
n is the exponent of the power-law fit below T c , and A e is the parameter in the exponential fit 
below T c that is proportional to the energy gap. The errors in the parameters were determined in 
the same manner as Tabled by propagating an error in Qd-, with the exception of T c and AC/T C 
(whose errors are represented graphically; see text for details). 
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FIG. 2: (a) Superconducting critical temperature T c vs Ru concentration x for Pr(Osi_ x Ru x )4Sbi2, 
including data from measurements of p(T), x ac (T), and C(T). The straight lines are guides to the 
eye. (b) The splitting between the ground state and first excited state AE gs -i cs vs Ru concentration 
x for Pr(Osi_ x Ru a: )4Sbi2, calculated from fits of CEF equations to Xdc{T), p(T), and C(T) data. 
The data points derived from p(T), Xa,c(T), and Xdc(^) are from previous workJ£ 
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FIG. 3: The electronic specific heat divided by temperature C e \/T (C c \ = C — C\ at — Cs c h) of 
Pr(Osi_ x Ru x ) 4 Sbi2 below 2 K for (a) < x < 0.2 and (b) 0.4 < x < 1. The data that were 
subtracted came mostly from fits including a Ti ground state. The data from subtracting a 1^ 
ground state fit were only included for x = due to large differences that are not present for other 
concentrations. 
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FIG. 4: Comparison of exponential (dashed line) and power law (solid line) fits, subtracting data 
for either 1^ ((a) and (c)) or T\ ((b) and (d)) ground states below T c , for PrOs4Sbi2 ((a) and (b)) 
and Pr(Oso.95Ruo.o5)4Sbi2 ((c) and (d)). The PrOs4Sbi2 fits only extend up to ~ 1.2 K, and the 
Pr(Oso.95Ruo.o5)4Sbi2 fits only extend up to ~ 1.1 K, as described in the text. 




FIG. 5: Electronic specific heat coefficient 7 (left axis, closed circles and squares) and AC/T C (right 
axis, open circles) as a function of Ru concentration x for Pr(Osi_ a; Ru a ;)4Sbi2. 
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